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O BLAZHKO EFFECT IN DOUBLE MODE CEPHEIDS 

O' 

CN ■ p. Moskalik^ and Z. Kolaczkowski^ 



__. , Abstract. Systematic survey for multiperiodicity in the LMC Cepheids 

(Moskalik, Kolaczkowski & Mizerski 2004! I2006|) has led to discovery 

of several new forms of pulsational behaviour. One of them is periodic 

1-^ . amplitude and phase modulation observed in many first/second over- 

Mh' tone (FO/SO) double mode Cepheids. In the current paper we present 

Q , detailed discussion of this newly discovered phenomenon, based on a 

Vh ' combined OGLE+MACHO sample of double mode pulsators. 

,^, 

1 Time Series Analysis 

J> ' We searched for additional signal in the data using a standard consecutive pre- 
l/^ , whitening technique. To that effect, we first fitted the data with the double fre- 
' ' ' quency Fourier sum representing pulsations in two radial modes: 

o 



o 

oo 

o 



m{t) = (to) + ^ k.jk sin[27r(jfi + fcfa)* + ^^tc]- (1) 



The frequencies of the modes, fi and £2, were also optimized. The residuals of the 
fit were then searched for additional periodicities. This was done with the Fourier 

^ ' transform, calculated over the range of — 5 c/d. In the next step, a new Fourier 
k>( fit with all frequencies discovered so far was performed and the fit residuals were 

; I ' searched for additional frequencies again. The process was stopped when no new 

Cu \ frequencies were detected. 

2 LIVIC Sample 

Our initial analysis was performed with DIA-reduced OGLE-II photometry (Ze- 
bruh et al. I^DUTj) . There are 56 FO/SO double mode Cepheids in the OGLE 
catalog (Soszyirski et al. I2000p . In one-third of these stars we detected residual 
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Table 1. Blazhko FO/SO Double-Mode Cepheids in the LMC 



^1 ^1 ^2 ^2 

OGLE ID/ Pi P2 Pb Ai A+ A++ A2 A+ A++ 

/MACHO ID [day] [day] [day] [mag] [mag] [mag] [mag] [mag] [mag] 

SCl-44845 0.9520 0.7660 794.0 0.165 0.021 — 0.046 0.021 — 

0.011 — 0.010 — 

SCl-285275 0.8566 0.6892 891.6 0.177 0.031 — 0.040 0.022 0.008 

0.014 — 0.013 — 

SCl-335559 0.7498 0.6036 779.2 0.227 0.019 — 0.062 0.023 — 

SC2-55596 0.9325 0.7514 768.2 0.146 0.007 — 0.039 — — 

0.007 — — — 

SC6-142093 0.8963 0.7221 1101.6 0.154 0.028 — 0.043 0.017 0.012 

0.013 — 0.017 0.006 

SC6-267410 0.8885 0.7168 856.9 0.120 — — 0.036 — — 

— — 0.033 — 

SC8-10158 0.6900 0.5557 1060.7 0.175 0.011 — 0.031 0.019 — 

— — 0.012 — 

SCll-233290 1.2175 0.9784 1006.2 0.186 0.019 — 0.043 0.017 0.006 

0.010 — 0.010 — 

SC15-16385 0.9904 0.7957 1123.1 0.258 0.017 — 0.049 0.020 — 

0.016 — — — 

SC20-112788 0.7377 0.5943 1379.2 0.164 0.062 0.022 0.015 0.016 — 

0.019 — 0.013 — 

SC20-138333 0.8598 0.6922 795.0 0.189 0.013 — 0.060 0.018 — 

0.008 — 0.010 — 

2.4909.67 1.0841 0.8700 1019.7 0.216 0.012 — 0.055 0.013 — 

0.010 — 0.009 — 

13.5835.55 0.8987 0.7228 1074.9 0.244 0.037 0.010 0.040 0.026 0.012 

0.021 — 0.016 — 

14.9585.48 0.9358 0.7528 1092.5 0.139 0.024 — 0.036 0.025 0.012 

0.008 — 0.010 — 

17.2463.49 0.7629 0.6140 1069.9 0.235 0.018 — 0.051 0.013 — 

0.014 — 0.015 — 

18.2239.43 1.3642 1.0933 706.8 0.230 0.035 0.010 0.063 0.017 — 

0.014 — 0.015 — 

22.5230.61 0.6331 0.5101 804.3 0.221 — — 0.045 0.016 — 

— — 0.014 — 

23.2934.45 0.7344 0.5918 797.6 0.207 0.023 — 0.055 0.031 0.015 

— — 0.026 0.011 

23.3184.74 0.8412 0.6778 1126.0 0.160 0.015 — 0.048 0.011 — 

— — 0.013 — 

80.7080.2618 0.7159 0.5780 920.3 0.168 0.010 — 0.055 0.006 — 

— — 0.009 — 
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Fig. 1. Power spectrum of FO/SO double mode Cepheid LMC/SCl-44845 after 
prewhitening with frequencies of the two radial modes (and their linear combinations). 
Removed radial frequencies are indicated by dashed lines. Lower panels display the fine 
structure around the radial modes. 



signal close to one or both primary (radial) pulsation frequencies. In all cases, how- 
ever, the secondary frequencies were not resolved from the primary ones. Clearly, 
~ 1200 days timebase provided by OGLE-II data was insufficient for our purpose. 
Therefore, the final frequency analysis was performed with MACHO photometry 
(AUsman & Axelrod 1200 ip . which had lower precision, but offered much longer 
timebase (2700 - 2800 days). 

The OGLE-II sample was supplemented by 51 additional FO/SO double mode 
Cepheids discovered by MACHO team outside the area covered by OGLE-II (Al- 
cock et al. IT9991 12003]) . Our final FO/SO Cepheid sample consisted of 107 objects. 

3 Results 



Residual power close to the primary pulsation frequencies was detected in 37 
FO/SO double mode Cepheids (35% of the sample). In 20 of these stars, which 
constitute 19% of the entire sample, we were able to resolve this power into indivi- 
dual frequencies. We consider two frequencies to be resolved if 1/Af < 1400 days. 
All resolved stars are listed in Table[T] 
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Fig. 2. Same as Fig.[T] but for LMC/SC6-142093. 



3.1 Frequency Domain 

Resolved FO/SO double mode Cepheids display very characteristic frequency pat- 
tern. In most cases, we detect two peaks in the vicinity of each radial frequency. 
They are located on the opposite sides of the primary (radial) peak and together 
with the primary frequency they form an equally spaced frequency triplet. In ma- 
jority of cases, such a structure is detected around both radial modes. Typical 
example of this pattern is displayed in Fig.[TJ In two Cepheids we detect equally 
spaced frequency quintuplets centered on the second overtone. One of these objects 
is presented in Fig. [21 Incomplete quintuplets are also seen in 6 other stars. 

The secondary peaks are always very small, with amplitudes of only 15mmag 
on average and never higher than 62 mmag. The separation of components in 
triplets/quintuplets, Af, is also very small and never exceed 0.00142 c/d. This 
corresponds to the beat period of Pb > 700 day. The lower limit for Af (upper 
limit for Pb) is currently given by resolution of the data. The physical upper limit 
for the beat period Pb is unknown. 

When present around both radial modes, the two triplets/quintuplets have 
nearly the same frequency spacings. Specifically, the difference between the two 
spacings, |Af2— Afi| is always below 9 x 10~^ c/d. According to simulations of Al- 
cock et al. (|2000l their Fig. 14), for MACHO data such a frequency difference is sta- 
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Fig. 3. Periodic moduiation of FO/SO double mode Cepiieid LMC/SCl-285275. First 
and second overtones are displayed with filled and open circles, respectively. Modulation 
period is Pb = 891.6 day. 



tistically consistent with no difference at all. In other words, triplets/quintuplets 
around both radial modes have the same frequency separation within accuracy of 
the data. 



3.2 Time Domain 

The triplet/quintuplet frequency pattern can be interpreted as a result of periodic 
modulation of amplitudes and/or phases of both radial modes, with common pe- 
riod Pb = 1/Af. In Fig. [3] we show this modulation for one of the stars. The plot 
has been constructed by dividing the data into 10 subsets, each spanning 10% of 
the modulation cycle, and then fitting radial modes and their linear combinations 
(Eq.(l)) to each subset separately. The resulting amplitudes and phases of the 
modes are plotted vs. modulation phase. 

Fig. [3] shows that both amplitudes and phases undergo modulation. The am- 
plitude variability is much stronger for the second than for the first overtone (80% 
vs. 30% change). Minimum amplitude of one mode coincides with maximum am- 
plitude of the other mode. Maximum phase of the first overtone occurs just before 
maximum of its amplitude. For the second overtone, maximum phase coincides 
with minimum amplitude. 
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OGLE LMC/SC17-186Q42 
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Fig. 4. Slow amplitude and phase variability of FO/SO double mode Cepheid 
LMC/SC17-186042. First and second overtones are displayed with filled and open circles, 
respectively. 

Characteristic pattern of periodic modulation displayed in Fig. [3] is common to 
all Cepheids listed in Table[TJ The phenomenon is similar to the Blazhko effect in 
RR Lyrae stars {e.g. Kurtz et al. I2000p . By analogy, variables of Table[l]can be 
called "Blazhko double mode Cepheids". 



3.3 Stars with Unresolved Power Spectrum 

In 17 FO/SO double mode Cepheids (16% of the sample) we find after prewhitening 
a significant residual power unresolved from the primary pulsation frequencies. 
This is a signature of slow amplitude and/or phase variability, not resolved within 
the length of available data {i.e. ~ 2700 days). 

For closer examination of this slow variability we have performed for each star a 
similar analysis to that presented in Fig.[3l Specifically, we have divided the data 
into 10 subsets, each spanning 10% of the total timebase, and then have fitted 
Eq.(l) separately to each subset. In some stars we have detected only changes of 
pulsation phases, but in 10 variables we have also detected clear changes of the 
amplitudes. In Fig.|4]we show typical example of such behaviour. The pattern of 
amplitude and phase variability is very much the same as that displayed in Fig. [31 
In particular, the amplitudes of the two radial modes vary in opposite directions. 
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Table 2. Unresolved FO/SO Double-Mode Ceplieids in the LMC 



MACHO Star 


Pi [day] 


Pa [day] 


OGLE Star 


A [day] 


P2 [day] 


6.6934.67 


0.9202 


0.7400 


SC4-176400 


1.1089 


0.8951 


11.9348.78 


0.7378 


0.5944 


SC4-220148 


0.7408 


0.5951 


14.9098.35 


0.7278 


0.5865 


SC7-120511 


1.2511 


1.0037 


15.10428.60 


0.6518 


0.5253 


SClO-204083 


0.5263 


0.4229 


23.3061.82 


0.5957 


0.4810 


SC16-266808 


1.3529 


1.0810 


24.2853.69 


0.7088 


0.5717 


SC17-186042 


0.6110 


0.4912 


24.2855.80 


0.5938 


0.4783 


SC20-188572 


1.0152 


0.8145 


47.2127.102 


0.5785 


0.4666 


SC21-12012 


1.3415 


1.0749 


80.7079.62 


1.3479 


1.0764 









This similarity suggests, that at least some of the unresolved FO/SO double mode 
Cepheids experience the same type of periodic modulation as Cepheids listed in 
Table[Tl but on timescales longer than the current data. 

4 What Causes the Modulation ? 



Periodic amplitude and phase modulation in FO/SO double mode Cepheids is by 
no mean a rare phenomenon. According to our analysis its incidence rate is at 
least 19% and perhaps as high as 35%. 

Any model of modulated double mode Cepheids has to explain two most basic 
properties of these stars: 

• both radial modes are modulated with the same period. 

• amplitudes of the two modes change in opposite phase: maximum of one 
amplitude coincides with minimum of the other. 

Two models have been proposed to explain Blazhko modulation in RR Lyrae stars: 
the oblique magnetic pulsator model (Shibahashi 11995] I2000p and 1:1 resonance 
model (Nowakowski & Dziembowski i200ip . We will argue, that both these models 
fail to account for modulation observed in the FO/SO double mode Cepheids. 

The oblique magnetic pulsator model assumes presence of a dipole magnetic 
field, which is inclined to the rotation axis of the star. The field introduces dis- 
tortion to the radial pulsation. As the star rotates, this distortion is viewed from 
different angles, causing variation of the pulsation amplitude. However, in this 
scenario both modes should display maximum amplitude simultaneously, because 
they are both distorted in the same way. This is not what we observe. 

The model proposed by Nowakowski & Dziembowski (|200ip assumes resonant 
coupling of a radial mode to nonradial modes of £ = 1. Such a mechanism genera- 
tes in a natural way a triplet of equally spaced frequencies. The beating of these 
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frequencies causes apparent amplitude and phase modulation. Generalizing this 
model to coupling with modes of £ = 2, we can also explain frequency quintuplets. 
However, in this resonance model, modulation of each radial mode is an indepen- 
dent effect. Therefore, there is no reason why both modes should be modulated 
with the same period. There is also no reason why amplitude variations of the 
modes should be in any particular phase relation to each other (as is observed). 
Finally, the 1:1 resonance model does not explain why in vast majority of cases 
modulation is observed either in both radial modes or in none of them. 

At this stage, the mechanism causing amplitude and phase modulation of 
FO/SO double mode Cepheids remains unknown. Nevertheless, the available ob- 
servations already provide some important clues. The fact that the two modes 
are always modulated with the same period proves that their behaviour is not 
independent. Both modes must be part of the same dynamical system. The fact 
that high amplitude of one mode always coincides with low amplitude of the other 
mode strongly suggests that energy transfer between the two modes is involved. 
Thus, available evidence clearly points towards some form of mode coupling in 
which both radial modes take part. 
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